Orbital Magnetic Susceptibility of Disordered Mesoscopic Systems 
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In this paper we study the orbital weak-field susceptibility of two-dimensional diffusive mesoscopic 
systems. For the previously unstudied regime of temperatures lower than the mean level spacing we 
find unexpected strong temperature as well as statistical ensemble dependence of the average and 
typical susceptibilities. An explanation for these features is given in terms of the long tail of the 
zero-temperature susceptibility distribution, including the parametric form of the temperature de- 
pendence. For temperatures higher than the mean level spacing we calculate the difference between 
the true canonical ensemble and the effective grand-canonical ensemble. We also perform numerical 
simulations, which seem to generally confirm previous theoretical predictions for this regime of tem- 
peratures, although some difficulties arise. The important role of gauge-invariance, especially how it 
renders Random Matrix Theory inapplicable to the study of orbital susceptibility, is discussed. We 
conclude by considering interaction effects, giving a new interpretation to previous results as well as 
demonstrating the influence of in-plane magnetic field on the interaction-induced orbital response. 

PACS numbers: 73.21.-b, 73.23.Ra, 75.75.+a 



I. INTRODUCTION 

Orbital magnetism is a purely quantum phenomenon, 
completely absent in classical physics, due to a theorem 
of Born and van-Leeuweni. R is therefore not surpris- 
ing that the magnetic response of mesoscopic systems, 
in which the particles move phase-coherently across the 
sample, shows many peculiar features, and is generally 
much larger than the magnetic response of macroscopic 
systems. 

Although this issue has been the subject of an intense 
scrutiny for almost two decades, several issues remain 
obscure. Since the problem was examined almost ex- 
clusively using diagrammatic perturbation expansion for 
diffusive systems and semiclassical analysis for ballistic 
ones, the behavior of the susceptibility at temperatures 
lower than the mean level spacing was not adequately 
understood. The consideration of the role of the differ- 
ent statistical ensembles is also incomplete, as we dis- 
cuss below, even for temperatures higher than the mean 
level spacing. Moreover, although much numerical effort 
was concentrated on investigating the magnetic response 
of multiply-connected geometries (persistent current), no 
similar effort was directed into understanding diffusive 
singly-connected geometries. Our aim in this paper is to 
give some theoretical considerations for the above men- 
tioned issues, and also to compare these and previous 
results with numerical simulations. 

We will start with a survey of previous results. Here, 
as well as in our analytical and numerical calculations, 
we limit ourselves to the case of singly-connected two 
dimensional systems of spinless noninteracting electrons 
(except in section^, and to the linear response regime. 
We will always consider the susceptibility per unit area. 

The orbital magnetic response of macroscopic systems 
is the celebrated Landau diamagnetic susceptibility 2 , 
X L = — e 2 /247TTOc 2 . This form is restricted to the linear 
response regime, Tiuj c < T, where T is the temperature 
(we use throughout units where the Boltzmann constant 



is unity) and lo c — eH/mc is the cyclotron frequency (H 
denotes the applied magnetic held) . The Landau suscep- 
tibility is temperature independent, up to temperatures 
of the order of the Fermi energy ep, above which the 
above result is multiplied by ep/T. 

This macroscopic Landau susceptibility is quite small. 
In a semiclassical picture this is a consequence of the ab- 
sence of flux-enclosing trajectories in infinite clean sys- 
tems. However, in mesoscopic systems such trajecto- 
ries do exist, due to scattering by either impurities or 
the sample boundaries. These two mechanisms tend to 
strongly enhance the susceptibility as described below. In 
the following we will assume that the dephasing caused 
by inelastic scattering (by other electrons, phonons, etc.) 
is negligible, i.e., > L, where is the phase coher- 
ence length, and L is the system's size. 

We will first consider disorder-induced susceptibility, 
which is dominant in the diffusive regime, kp 1 < I < L, 
where kp is the Fermi wave- vector and I is the elas- 
tic mean free path. The presence of disorder causes 
large sample specific fluctuations. Thus, the typical sus- 
ceptibility, measured by the r.m.s. value 8x, is much 
larger than the average susceptibility, where the aver- 
age is over an ensemble of disorder realizations. For 
temperatures below the Thouless energy Ec = TiD/L 2 
(D = vp£/2 is the diffusion coefficient in two dimen- 
sions), the typical susceptibility is expected to vary 
as |x L |(M) [H E c/T)} 1/2 . The form of the logarith- 
mic factor is correct in the very low field regime $ < 
(T/Ec^^^o, where <f> = HS is the total magnetic 
flux through the system (S is the system's area), and 
<3>o = hc/e is the flux quantum^. These fluctuations 
may persist even in temperatures above the Thouless en- 
ergy, where the temperature dependent factor changes to 
(Ec /T) 1 / 2 , with an increased range of validity for linear 
response^. 

We now regard the disorder averaged susceptibility. 
Here we should distinguish between the Canonical En- 
semble (CE) and the Grand Canonical Ensemble (GCE). 
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The GCE corresponds to a situation where the systems 
in the ensemble of disorder realizations can exchange par- 
ticles with a heat bath with a fixed, system independent 
chemical potential. The CE corresponds to a situation 
where the systems are isolated, with the same fixed num- 
ber of particles in each. Usually it was assumed that the 
CE can be described by an Equivalent Grand Canonical 
Ensemble (EGCE), in which each system can exchange 
particles with a different heat bath, with its chemical po- 
tential adjusted to the realization so that the thermally 
averaged number of particles in the system equals the 
given fixed particle number. Thus, the difference appears 
only after averaging over an ensemble of systems with dif- 
ferent disorder realizations, but not for a single system. 
However, it was argued several times in the past that the 
CE and EGCE are really equivalent only at zero tem- 
perature or at high temperatures (higher than the mean 
level spacing A or even the Thouless energy Ec)- For in- 
termediate temperatures the level occupancy in the CE 
cannot be described by the Fermi-Dirac distributioni£*il. 

It is predicted that in the GCE the ensemble aver- 
aged susceptibility will be equal to the Landau diamag- 
nctic response. Disorder only modifies this result by 
a small contribution, of order (kp£)~ . On the other 
hand, in the EGCE there is an additional paramag- 
netic contribution^, of order |x L |(fcp£)A/T. This con- 
tribution is much larger than the Landau susceptibil- 
ity, and will dominate both the sign and the magni- 
tude of the response up to temperatures of the order 
of the Thouless energy Ec- This paramagnetic part of 
the response is expected to be strongly field dependent 
even in the weak field regime $ < $o , and linear re- 
sponse is expected to apply only for very low fields, for 

which $ < (T/E c ) 1/2 <Z>o- As we will show in the ap- 
pendix, moving to the pure CE adds another term, of 
order |x L \(kp£)(A/T) 2 , which is smaller than the differ- 
ence between the GCE and the EGCE, but may still be 
larger than the Landau susceptibility. All this applies 
to the average susceptibility; the r.m.s. susceptibility is 
always assumed to be ensemble-independent, at least to 
the leading order. 

In the past few years it was recognized that the 
disorder-induced susceptibility can be important even in 
the ballistic regime £ > L, and will cease to exists only 
in the clean limit, £ > (kpL)L, where the disorder causes 
only a small perturbation to the energy levels. However, 
the situation in this case is still unclea»i2ii£, and will not 
be discussed in the following. 

It should be noted that all the quoted results in the dif- 
fusive regime were obtained using diagrammatic pertur- 
bation theory, and are therefore valid only for tempera- 
tures above the mean level spacing A. It was conjectured, 
using different semi-quantitative arguments, that there is 
almost no temperature dependence for T < A, and that 
in this regime the values of the average and typical sus- 
ceptibility are given by the above expressions evaluated 
at T ~ A^. 

In clean mesoscopic systems the susceptibility is en- 



hanced by boundary effect s^Si Here, the spatially lim- 
ited electron motion results in oscillations in the suscep- 
tibility as a function of the electron density. The oscilla- 
tion's period is of order kp L, where kp is the Fermi wave- 
vector. These oscillations are accompanied by large para- 
magnetic peaks as a function of electron density. These 
peaks result from the occurrence of degenerate or nearly- 
degenerate levels, and are eliminated as the temperature 
T increases above the mean level spacing A. The oscilla- 
tions, however, persist up to much higher temperatures, 
of the order of the inverse time of flight across the sam- 
ple, Hvf/L, beyond which they decay exponentially with 
the temperature (vp is the Fermi velocity). In all the 
cases the susceptibility is much larger than the Landau 
value, and is of order |x L |(fcpL) a , where the exponent 
a is of order unity, depending on whether the system is 
classically integrable or chaotic, and also on the specific 
geometry in the former case. Introducing disorder into 
the system makes these effects decay as exp(—L/£), i.e., 
they should disappear in the diffusive regime £ < L. 

Since in reality the system's size and shape cannot be 
set accurately, one usually averages the results over a 
distribution of sizes around an average value. This size 
averaging smears out the oscillations, and reduces the 
size averaged GCE susceptibility to the Landau value for 
temperatures higher than the mean level spacing. The 
size averaged EGCE susceptibility, as well as the size 
averaged r.m.s. susceptibility remain enhanced by factors 
which are again powers of kpLm£ In this paper we will 
use the clean case only for comparison with the diffusive 
case, so we will not perform size-averaging. 

The paper is organized as follows: In section[H]we will 
first show that the source of difficulties in evaluating the 
response of diffusive systems in the range T < A is re- 
lated to the problem of keeping the gauge-invariance of 
the results in the disorder averaging, and that this makes 
Random Matrix Theory (RMT) inapplicable for evaluat- 
ing the susceptibility. This discussion will also serve to 
introduce notations and formulas for usage in subsequent 
calculations. In section ITTT1 we will then examine our nu- 
merical results for both clean and diffusive systems in the 
various statistical ensembles. These will be seen to gen- 
erally agree with the theoretical predictions for T > A, 
although not in all the details. For T < A an unexpected 
strong temperature and statistical ensemble dependence 
of the typical susceptibility in the diffusive regime is ob- 
served, together with similar but weaker effects in the 
average susceptibility. This will be shown in section llVl 
to stem from long tails in the susceptibility distribution. 
It will enable us to work out the parametric form of the 
r.m.s. susceptibility temperature dependence, as well as 
the reasons for its statistical ensemble dependence. Fi- 
nally, we will consider interaction effects in section Ivl 
We will give a new way of interpreting the reasons for 
interaction-induced susceptibility, as well as a numeri- 
cal demonstration of its dependence on applied in-plane 
magnetic field. 
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II. CONSEQUENCES OF GAUGE INVARIANCE 

We consider a system of spinless non-interacting elec- 
trons confined to a finite area and moving in a random 
potential of impurities and a magnetic field. The system 
is described by the following Hamiltonian: 

H=^- i (p- e -A) 2 + V conl + V disi (1) 

where p is the momentum operator, A is the vector- 
potential, y c onf is the confining potential and Vdis is the 
random impurity potential. 

If the magnetic field is weak, its influence can be taken 
as a perturbation. Since the eigenstates in zero field can 
be chosen real and are non-degenerate (due to level re- 
pulsion), the first order correction to the energy (contain- 
ing the term linear in the magnetic field in the Hamil- 
tonian to first order in perturbation theory) identically 
vanish, i.e., there is no spontaneous magnetic moment. 
The second order correction to the energy determines 
the susceptibility, which is composed of the familiar Lar- 
mor diamagnetic term (containing the term quadratic in 
the magnetic field in the Hamiltonian to first order in 
perturbation theory) and van-Vleck paramagnetic term 
(containing the term linear in the magnetic field in the 
Hamiltonian to second order in perturbation theory) . At 
finite temperature they take the following form : 

*- = -5^£/(-)<HA 2 | fc >, 

k 

v _ 1 ( e Y v /( efc )-/( £i ) 
XvV SH 2 \2mcJ £-< e k -ei 

x|(fc|A-p + p.A|0| 2 , (3) 



where the indices a, (3 run over the two Cartesian coordi- 
nates and the integrations are over the system's area. 
Since the wave functions in this expression are most 
strongly correlated when r = r', but the integrand van- 
ishes there identically, we can understand why the re- 
sulting average susceptibility is much smaller than the 
Larmor or van-Vleck terms separately. 

Another lesson we can learn from this is that when 
trying to calculate the disorder averaged susceptibility, 
one should check this does not destroy the gauge invari- 
ance of the results. For example, the proof of gauge in- 
variance in our case relies on the relation (k\p/m\l) = 



where fc) , \l) are the eigenstates of the Hamiltonian 
without magnetic field and /(e) is the average occupa- 
tion of a level with energy e, whose form depends on the 
required ensemble, and will be given at the end of this 
section. 



It has been noted several times in the past^ii*i^ that 
there is a large cancellation between these two contri- 
butions. We interpret this as a consequence of gauge 
invariance. If one replaces the vector potential A by 
A + V</>, where <fi is some function of the coordinates, the 
observables of the system cannot change. When treating 
the magnetic field as a perturbation, this means that <p 
must cancel between all the terms of the same order in 
the perturbation theory. Indeed, it can be shown that 
the changes of the Larmor and van-Vleck susceptibili- 
ties due to the gauge transformation identically cancel. 
This explains why these two seemingly separate contri- 
butions are strongly correlated and tend to cancel each 
other, leaving behind only the gauge invariant part of 
their sum. 



In fact, the expression for the Larmor susceptibility 
can be brought into a form very similar to the expression 
for the van-Vleck susceptibility, thus demonstrating the 
deep connection between them. The resulting formula 
for the total susceptibility reads (in the family of gauges 
obeying V • A = 0) : 



(4) 



i(ek — £i)/h{k |f 1 1). This implies a connection between 
matrix elements of the eigenstates and their eigenvalue 
difference which is usually absent in many averaging 
schemes, especially RMT. 

This problem does not arise in diagrammatic calcu- 
lations, where the average or r.m.s susceptibilities are 
expressed in terms of the energy shifts of the diffuson 
and Cooperon propagators, since the equations govern- 
ing them are manifestly gauge invariant. Neither is there 
a problem in the semiclassical approach, where the influ- 
ence of weak magnetic fields is taken only through the 
phase accumulated in closed orbits, which depends only 



2SH 2 



2 

e" 



h\ ^ m - /fc) f dv f dr ^ (A a (r) - A a (r'))(A,(r) - A,(r')) 

10 J k?l £k ~ ei Js Js a,0 

xV fe (r)V Q 7/; ^ (r)V' ^ (r')V , /3 V'fc(r , ), 
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on the total flux through the orbits and is thus again 
gauge invariant. However, for temperatures lower than 
the mean level spacing those approaches are inapplicable. 
In this regime RMT is usually used, but, as mentioned 
above, it fails to give gauge invariant results. 

To conclude this section we give the explicit form of the 
mean level occupation f(e n ) of the n-th single particle 
level. In the GCE or EGCE it is simply the Fermi-Dirac 
distribution function, 



/gce(c) = 



.(e-/*)/T + !- 



(5) 



where the chemical potential \i is system independent in 
the GCE. In the EGCE it is system specific, determined 
by the requirement that the thermally averaged particle 
number in the EGCE, J2 n f( e n)i will be equal to the con- 
stant particle number in the sample. For the CE the oc- 
cupation function can be represented by Darwin- Fowler 
integrals, which are convenient for both analytical and 
numerical calculations: 



/ce(c) 



ccEie; z e 



-Vt(z)/T dz 



,N+1 



f 



-Q(z)/T dz 



(6) 



where 



Q(z) = -T^lnfl 



is the grand canonical potential, z = exp(/x/T) is the fu- 
gacity, and the integrations are taken around any closed 
contour in the complex z plane encircling the origin. This 
expression can be shown to reduce to the Fermi-Dirac 
distribution function both for zero temperature and for 
temperatures much higher than the mean level spacing, 
where the integral can be evaluated using saddle-point 
approximation (The saddle-point is the fugacity corre- 
sponding to the sample-specific chemical potential in the 
EGCE). 



III. NUMERICAL SIMULATIONS 

In this section we will present the results of our numer- 
ical simulations for the susceptibility. For this purpose 
we use the familiar Anderson tight-binding Hamiltonian: 



H 



E 



t 



J6 



<s,s'> 



'a\a s 



(7) 



Anderson model has the advantage of being discrete, but, 
unlike numerical discretization of a continuous Hamilto- 
nian, it has all the required properties of a Hamiltonian 
(i.e., being Hermitian, gauge- invariant, etc.) 

In the calculation we use expressions for the Larmor 
and van-Vleck susceptibilities for the Anderson Hamil- 
tonian. Since the magnetic field is assumed small, we 
first expand the phase exponent containing the magnetic 
field in Eq. Q in series. As for the continuous Hamilto- 
nian Q , the linear term in the field taken to first order in 
perturbation theory identically vanishes. We are left with 
the second order corrections to the energy - the quadratic 
term in the field taken to first order in perturbation the- 
ory (Larmor susceptibility), and the linear term in the 
field taken to second order in perturbation theory (van- 
Vleck susceptibility) .We thus get formulas analogous to 
Eqs. ©,© f° r the susceptibility, which are used in all 
the subsequent calculations. The level occupations are 
calculated in the required ensembles, using Eqs. 0,|JBJ. 
We note that even though the contour of integrations in 
Eq. 10 can be evaluated in principle using any path in 
the complex plane encircling the origin, to get sensible 
results in numerical integration one should use a contour 
passing through or near the sample specific saddle point 
(the EGCE fugacity) in a direction where this point is a 
maximum of the integrand. A suitable form of the con- 
tour is thus a circle around the origin, and this choice 
was used in our calculations. 

The calculations were usually made (unless otherwise 
specified) on a 17 x 24 lattice (the sizes were chosen mu- 
tually prime so that there's no degeneracy even in the 
clean limit, and our non-degenerate perturbation theory 
is applicable). We have made calculations both on clean 
(W — 0) systems and disordered systems in the diffu- 
sive regime, especially with the disorder values W = 2.0 
and W — 4.0. For the diffusive systems the results were 
averaged over 2500 realizations of disorder, unless oth- 
erwise specified. Dependence on the system's size was 
deduced from comparison to the results on 13 x 19 and 
8 x 13 lattices. The susceptibility is usually plotted as a 
function of electron filling (i.e., the ratio of the number 
of electrons and the number of lattice sites). To enable 
comparison this is done not only in the CE or EGCE, but 
also in the GCE, in which case we refer to the ensemble 
and thermally averaged filling. 



A. Macroscopic Clean Systems 



where a\, d s and n s denote electron creation, annihila- 
tion and number operators, respectively, for a state on 
site s of a square lattice. The first term is a random on- 
site potential, where e s is chosen randomly in the range 
[-W/2, W/2]; the second is the hopping or kinetic term, 
where the sum is over nearest-neighbor sites s and s', t is 
the overlap integral, and the phase 9 S , S > = e/Kc f s , A • dr 
gives the influence of the external magnetic field. The 



For comparison with subsequent results for finite sys- 
tems, we first show the orbital susceptibility of the tight- 
binding Hamiltonian J7J for infinite clean lattice. The 
zero temperature susceptibility was calculated in Ref. 0, 
and is given by the following expression : 



x(T = 0,e F ) 



2Ah 2 c 2 



[(e 2 F ~8t 2 ) p(e F )-E(e F )] ,(8) 
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FIG. 1: The orbital susceptibility of an infinite ordered lattice 
as a function of filling for various temperatures. 



where 



p(e) 



2ir 2 t 



K 



\Jl-(e/4tf 



is the density of states (K denotes the complete elliptic 
integral of the second kind), E{ep) is the total energy 
of the filled states, and a is the lattice constant. At the 
band edges we get the value —a 2 e 2 t/12irh 2 c 2 , which is 
the Landau susceptibility with the appropriate effective 
mass. At nonzero temperatures the susceptibility can 
be evaluated from its zero temperature value using the 
connection x(T,n) = J x(0, £)/gce( t > e ) de > where f^ CE 
is the energy derivative of the Fermi-Dirac distribution 
function (of course, in this macroscopic limit the differ- 
ence between statistical ensembles disappears). 

The susceptibility is plotted in Fig. ^ as a function of 
the filling v = J p(e)focE(e)de. We can see that the sus- 
ceptibility has a diamagnetic Landau value only at the 
band edges at zero temperature, and that its sign changes 
to paramagnetic near the band center. At half-filling it 
actually exhibits a logarithmic singularity at zero tem- 
perature. Those non-analytic behaviors at the band's 
center and edges cause noticeable temperature depen- 
dence of the susceptibility there, even for temperatures 
much lower than the Fermi energy (the only energy scale 
in macroscopic clean systems), although in the rest of the 
band temperature has a significant influence only when 
it reaches the Fermi energy. This behavior is connected 
to the inclusion of the magnetic field only as a phase fac- 
tor in the hopping amplitude. However, if we remember 
in the sequel to compare our results for finite, possibly 



1400- 
1200- 
1000- 



hJ 800- 



- T/t=0.0 
T/t =0.001 

- T/t =0.01 



600- 
400- 
200- 



0- 



i I , I i c jL 

0.1 0.2 0.3 

Electron Filling 



0.4 



0.5 



2 
1.5 
1 

0.5 



-0.5 



- T/t =0.1 
T/t=1.0 



' i i.-i r 



'\ 1 



"A '\' 'i 
' i ' V 1 1 



,-t-T 

\ i 
W 



0.1 0.2 0.3 0.4 

Electron Filling 



0.5 



FIG. 2: The orbital susceptibility of a finite clean lattice as a 
function of filling for various temperatures in the EGCE. Low 
temperature results are presented in the upper panel, high 
temperature results are presented in the lower panel. 



disordered systems to this macroscopic limit, we would 
not have any trouble interpreting the results. 



B. Mesoscopic Clean Systems 

We now turn to finite clean systems. The results for 
the susceptibility of a system on a 17 x 24 lattice are 
shown in Fig. in the EGCE (which is equivalent to 
the GCE in the absence of disorder), for various tem- 
peratures. The mean level spacing A is approximately 
0.02i. We can see here all the expected phenomena de- 
scribed in the introduction. At zero temperature, or 
finite temperatures much smaller than the mean level 
spacing, there are very large paramagnetic peaks, occur- 
ring where the level spacing is small. Their exact po- 
sition and height strongly depend on the system's area 
and the ratio of its length and width. Those peaks are 
smeared out when the temperature approaches the level 
spacing, leaving behind oscillations whose amplitude and 
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FIG. 3: The orbital susceptibility of a finite ordered lattice as 
a function of filling for various temperatures in the CE. Low 
temperature results are presented in the upper panel, high 
temperature results are presented in the lower panel. 
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FIG. 4: The disorder averaged orbital susceptibility of a fi- 
nite disordered lattice as a function of filling for various tem- 
peratures in the GCE, for both W = 2.0 (upper panel) and 
W — 4.0 (lower panel). 



frequency scale roughly as L 3 ' 2 and L~ x , respectively, as 
in continuous square geometries 6 . They still have a para- 
magnetic tendency, although much smaller than the low 
temperature peaks. These oscillations are smallest near 
the band edges or center, and are maximal near quarter 
filling. This is in contrast with the macroscopic suscep- 
tibility, which is maximal near the band center. The 
oscillations disappear, and the susceptibility assumes its 
macroscopic value (including the subsequent temperature 
dependence) for temperature around 0.2i w 10A, which 
can be attributed to the inverse time of flight across the 
sample hvp/L for all but the lowest fillings. 

As regards the CE, The results for the same values of 
parameters are shown in Fig. EI We can see that moving 
from the EGCE to the CE is qualitatively equivalent to 
lowering the temperature in the EGCE, as noted by var- 
ious authors in the pasfei^iii. Hence, the peaks and oscil- 
lations are somewhat larger in the CE, and thus persist to 
higher temperatures. The difference disappears at both 
zero temperature and very high temperatures (of the or- 



der of the temperature where the mesoscopic effects in 
the EGCE disappear), but is observable for intermediate 
temperatures. However, even when non-negligible, this 
difference is smaller than the EGCE susceptibility. 



C. Mesoscopic Disordered Systems, GCE 

We now introduce disorder into the system. The re- 
sults for the average susceptibility in the GCE are shown 
in Fig. The GCE susceptibility is expected, as noted 
in the introduction, to be equal to the macroscopic sus- 
ceptibility, shown in Fig.^ Indeed, the general behavior 
of the GCE susceptibility resembles the macroscopic sus- 
ceptibility (especially in its sign as a function of filling). 
This is immediately obvious regarding the case W = 4.0, 
and is also true for W = 2.0 if the oscillations as a func- 
tion of filling are averaged out (These oscillations will 
be discussed below). In addition, the GCE susceptibility 



(ignoring again the oscillating part for W = 2.0) shows 
no temperature dependence, up to about T « O.lt. This 
is in accordance with the macroscopic behavior if we re- 
member that there are no zero temperature singularities 
as a function of filling for finite samples, so only the high 
temperature macroscopic temperature dependence per- 
sists. 

There are, however, two modifications compared to the 
macroscopic results. The first is that the susceptibility of 
the disordered mesoscopic systems is clearly smaller than 
the macroscopic value. It also decreases when disorder 
increases, although quite slowly. This can be interpreted 
as stemming from corrections to susceptibility (beyond 
the leading order macroscopic value), which are expected 
to vary as (kp£)~ . Since (kp£)~ is not so small in 
our simulations, these corrections can be important and 
lead to a weak disorder dependence, and their sign is 
apparently negative. 

A second and more pronounced difference between the 
GCE average susceptibility and its macroscopic value 
is the above mentioned filling-dependent oscillations for 
W = 2.0. The oscillations' position is the same as the 
oscillations in the ordered systems shown above in Fig. |2 
i.e., of periodicity varying as L _1 . The amplitude of the 
oscillations, however, shows little dependence on the sys- 
tem's size, but decreases strongly when disorder or tem- 
perature increase. This effect is due to the fact that even 
though for W — 2.0 the systems of the size considered 
are regarded diffusive, showing, for example, a Wigner 
level spacing distribution, the elastic mean free path I 
is only about a third of the size of the considered sys- 
tems, so that ballistic effects can still be observed. In 
fact, the density of states still shows ballistic oscillations 
for W = 2.0, which are smeared out as disorder increases 
toward W = 4.0. Thus, the susceptibility at W — 2.0 
shows a mixture of diffusive and ballistic behaviors. We 
note that size-averaging will smear out these oscillations. 

Referring now to fluctuations in the susceptibility, the 
results for the r.m.s. susceptibility of the same system 
are shown in Fig. We can clearly see that the r.m.s. 
susceptibility is much larger than the average. It also 
decreases strongly with increasing disorder, implying the 
expected IXlI(^f^) dependence (This can be seen quan- 
titatively by noting that from Fermi's golden rule one 
can estimate I to vary as W~ 2 , which is quite accurately 
the W dependence of the r.m.s susceptibility in our data). 
Since kpl ~ g, where g is the dimensionless conductance, 
the enhancement factor agrees well with the values of 
g as 20 and g « 3 for W = 2.0 and W = 4.0 respectively. 
(The values g were calculated from non-universal correc- 
tions to the inverse participation ratio2I, and agrees well 
with estimates based on Fermi's golden rule.) 

A surprising feature of the results is that the r.m.s. sus- 
ceptibility depends strongly on temperature even for tem- 
peratures much smaller from the mean level spacing A 
(which is, approximately, 0.022t and 0.025t for W = 2.0 
and W = 4.0, respectively) This will be discussed later 
on. Moreover, the r.m.s. susceptibility becomes equal 
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FIG. 5: The disorder averaged r.m.s. orbital susceptibility 
of a finite disordered lattice as a function of filling for various 
temperatures in the GCE, for both W = 2.0 (upper panel) 
and W = 4.0 (lower panel). 



to the average susceptibility for a temperature of about 
O.lt « 5 A for both W values. As mentioned in the in- 
troduction, significant fluctuations in the susceptibility 
(i.e., larger than the average value) are expected to per- 
sist even in temperatures higher than the Thouless en- 
ergy. However, since the enhancement factor kF& is not 
so large in our simulations as in real systems, absence of 
fluctuations in temperatures above the Thouless energy is 
quite reasonable. On the other hand, fluctuations below 
the Thouless energy should be significant. This implies 
interpreting the temperature O.lt ~ 5 A, for which the 
difference between the average and r.m.s. susceptibilities 
disappears, as the Thouless energy times a numerical fac- 
tor (note that our estimate for g indicates a much higher 
Thouless energy for W = 2.0.) 

However, although this temperature goes as L~ 2 as 
the Thouless energy should vary, it shows very little if 
any disorder dependence, so its identification with the 
Thouless energy is quite problematic. This may support 
the claim of some authors^ that in the logarithmic factor 
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susceptibility equals the Landau value. It may be ex- 
plained by the fact that those states are more localized 
than the other states, and thus their susceptibility is not 
so much enhanced. 

The difference between the average EGCE and GCE 
susceptibilities is strongly disorder dependent, confirm- 
ing the identification of the zero temperature enhance- 
ment factor as fc^ (As for the r.m.s. GCE susceptibil- 
ity discussed above, it has quite accurately a W~ 2 de- 
pendence). Temperature dependence shows up for tem- 
peratures higher than for the r.m.s. GCE susceptibility 
discussed above, but much lower than the mean level 
spacing. In fact the difference between the EGCE and 
GCE practically disappears for T = 0.02< w A and 
T = O.Oli « 0.4A for W = 2.0 and W = 4.0 respec- 
tively. Although this temperature is much smaller than 
the Thouless energy (based on our estimate g w 20 and 
g w 3 for W — 2.0 and W — 4.0 respectively), its L~ 2 size 
dependence as well as its disorder dependence indicate 
that it might be considered as connected with the Thou- 
less energy reduced by a numerical factor. We may note, 
however, that the disorder dependence of this tempera- 
ture is quite weaker than the corresponding dependence 
of the Thouless energy (the above mentioned W~ 2 law). 
Combined with the quite low value of this temperature, 
its meaning is not completely clear. 

As for the fluctuations in the susceptibility, The re- 
sults for the r.m.s. susceptibility are shown in Fig. Q As 
is usually conjectured, it behaves quite similarly to the 
GCE r.m.s. susceptibility, and all the discussion regard- 
ing the GCE applies also here. The only exception is 
the very low temperature regime, (below O.Olt w 0.5A), 
where the EGCE r.m.s. susceptibility is well above the 
GCE values, the difference being larger for W — 2.0. It 
is also much more "noisy" , as can be clearly seen in the 
T = curves. We will discuss this feature in the next 
section. 



h\{Ec /T) in the r.m.s susceptibility, the Thouless energy 
should be replaced by some lower energy scale. Our re- 
sults seem to suggest that this scale, which is disorder 
independent, is connected to the mean level spacing A. 



D. Mesoscopic Disordered Systems, EGCE 

We now turn to discuss disordered systems in the 
EGCE. The results for the average susceptibility are 
shown in Fig. We see again the oscillations in the 
W = 2.0 susceptibility. Their explanation is the same as 
in the GCE case discussed above, so we now ignore them 
and refer to the smooth part of the susceptibility. 

As is expected, the difference between the average 
EGCE susceptibility and the GCE susceptibility is posi- 
tive. It is also much larger than the GCE susceptibility 
except near the band edges. This is somewhat problem- 
atic since precisely at the band edges the macroscopic 



E. Mesoscopic Disordered Systems, CE 

To conclude this section, we move to the CE. The re- 
sults for both the averaged and r.m.s. susceptibilities 
are very similar to the EGCE results, the difference be- 
ing only few percents of the EGCE values and thus not 
shown. This is quite expectable for T > A, since, as 
we show in the appendix, the difference between the CE 
and EGCE is parametrically smaller than the difference 
between the EGCE and the GCE. 

We may note, however, that the difference between the 
CE and EGCE is more pronounced for the r.m.s suscep- 
tibility at very low temperatures (below O.Oli 0.5A), 
where there is also a difference between the EGCE and 
the GCE r.m.s. susceptibilities. We believe those two 
phenomena has a similar origin, and will discuss it in the 
next section. 
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FIG. 7: The disorder averaged r.m.s. orbital susceptibility 
of a finite disordered lattice as a function of filling for various 
temperatures in the EGCE, for both W = 2.0 (upper panel) 
and W = 4.0 (lower panel). 



IV. ZERO TEMPERATURE SUSCEPTIBILITY 
DISTRIBUTION 

In this section we will examine the zero tempera- 
ture susceptibility distribution, and use it to explain the 
strong dependence of the typical and average suscepti- 
bilities on both temperature and the statistical ensemble 
for temperatures lower than the mean level spacing. 

The zero temperature susceptibility distribution in our 
systems is shown in Fig. EI f° r W = 2.0 and W = 4.0 and 
for both the CE or EGCE and the GCE. In all the cases 
the results are for half-filling. The numerically obtained 
distributions are shown together with Gaussian fits. The 
tails of the distributions are shown in the insets, together 
with fits to functional forms which we derive below. We 
can see that all the distributions are peaked at (or near) 
zero susceptibility, and are somewhat asymmetric around 
their maximum. They are, however, fitted quite well by 
a Gaussian near the peaks. On the other hand, The 
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FIG. 8: The orbital susceptibility distribution of a finite 
disordered lattice at zero temperature, for both the CE or 
EGCE and the GCE for half-filling, with disorder values W = 
2.0 (upper panel) and W = 4.0 (lower panel). The tails of 
the distributions are enlarged in the insets. Numerical results 
(histograms) are shown together with fits (continuous curves) . 
See the text for further details. 



tails are clearly non-Gaussian, decaying somewhat faster 
in the negative direction and much more slowly in the 
positive direction. We note that the susceptibility distri- 
butions for other fillings, where the average susceptibility 
is different, are similar in form to Fig. [SI but show the 
required tendency toward higher or lower values. 
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We now compare the differences between the severai 
cases shown. We can see that the CE distribution is wider 
than the GCE distribution for both disorder values, and 
has a stronger tendency for positive values. In particular, 
the positive tail of the CE distribution is larger than the 
corresponding tail in the GCE. When comparing the two 
disorder values, note we have chosen a factor of 4 between 
the scales of the two graphs. This factor is both as the 
ratio of the widths of the Gaussian fits and a consequence 
of the expected W~ 2 dependence of the susceptibility 
moments. We can see that after this rescaling the results 
for W = 2.0 are similar to those for W — 4.0, although 
the distributions are somewhat wider in the former case. 

In two recent papersii*i£ Serota has attempted ad- 
dressing the problem finding the zero temperature av- 
erage CE susceptibility and its fluctuations for diffusive 
systems. His discussion was based on the expressions for 
the Larmor and van-Vleck susceptibilities, Eq. and 
Eq. respectively. In order to take into account the 
above mentioned strong cancellation between the Larmor 
and van-Vleck terms, he conjectured that the only contri- 
bution surviving this cancellation is the largest van-Vleck 
term composed of a matrix element and an energy differ- 
ence between the highest occupied level and the lowest 
unoccupied level, which we will denote by — and + re- 
spectively. At zero temperature one is thus left with the 
following expression for the zero temperature susceptibil- 
ity: 

x - ' ( ej W.p + p.AI-if , (9) 
bti l \2vncJ e + — e_ 

Using an Gor'kov-Eliashberg 18 like argument to eval- 
uate the squared matrix element semiclassically, and 
using the Wigner distribution for the level spacing, 
Serota has given expressions for the average and 
r.m.s. susceptibilities, which vary as IXlK^f-O arid 

1/2 

|x L \(kp() \hi(Ec/A)] . These expressions are similar 
to the corresponding expressions for temperatures higher 
than the mean level spacing A which we quoted in the 
introduction, if the temperature T is replaced by A. 

However, this approach is quite problematic in view of 
our previous discussion of the importance of the gauge 
invariance of the susceptibility. Eq. J3J is clearly not 
gauge invariant, since the Larmor- van-Vleck cancellation 
was taken into account a is too crude fashion. As a re- 
sult, Eq. ©, can only give positive susceptibility values, 
whereas our numerically computed distribution, Fig. [S] 
exhibits also negative values. 

In spite of this problem, the approximation of Eq. [§] is 
still useful for treating the tail of the susceptibility distri- 
bution. Our numerical simulations indicate that anoma- 
lously large positive susceptibility values are caused by 
small level separation between the highest occupied level 
and the lowest unoccupied level. In those cases, the term 
in the van-Vleck susceptibility involving these two levels 
is much larger than all the other terms. Furthermore, in 
this case it is also approximately gauge invariant, since 



a gauge transformation adds to the squared matrix el- 
ement in Eq. © contributions smaller by a factor of 
(e + — e_)/A or less, which are thus negligible for very 
close levels. Moreover, the average and r.m.s. of the 
squared matrix element are seen in the numerics to be 
almost independent of the energy separation for small 
separations, so that the form of the tail of the suscep- 
tibility distribution is dominated by the distribution of 
the energy denominator. This, in turn, can be simply 
obtained from Wigner's distribution, according to which 
the probability for small level separation varies linearly 
with the separation. Thus, the susceptibility distribution 
varies as P(x)dx ~ x _3 dx- This dependence is fitted to 
the CE distributions in the insets of Fig. [8] and agrees 
quite well with the numerical results. 

As previously noted by Serota, this distribution has a 
first moment but not a second or higher moments, i.e., 
the r.m.s. susceptibility diverges. Serota suggested regu- 
larizing this divergence by noting that in any given exper- 
imental realization a finite magnetic field is used, so that 
for very small energy separations the non-degenerate per- 
turbation theory used to obtain Eq. © is inapplicable. 
A degenerate perturbation theory must be used instead, 
and this gives a finite magnetization. 

All this is true, however, at zero temperature. At 
higher temperatures Eq. @ should be multiplied by 
the factor (/(e_) — /(e+))/2, where /(e) is the mean 
level occupation. This factor tends to linearly zero for 
level separation smaller than the temperature, and thus 
cuts off the divergence caused by the denominator. For 
temperatures higher than hu c this factor, and not a fi- 
nite magnetic field, is dominant in keeping the suscep- 
tibility finite. In fact, using Serota's estimate for the 
squared matrix element in Eq. @, the r.m.s. suscep- 

1/2 

tibility can be shown to vary as IXlI(^f^) [ m (A/T)] • 
This means that the logarithmic temperature dependence 
of the r.m.s. susceptibility 8 does not saturate for tem- 
peratures smaller than the level spacing, but persists to 
much lower temperatures, and, in the limit of zero mag- 
netic field, continues down to zero temperature, where 
the r.m.s. susceptibility diverges. 

This may serve to explain the strong temperature de- 
pendence of the r.m.s. susceptibility seen in our numeri- 
cal results. It also suggests that the saturation of this ef- 
fect at very low temperatures (less than tenth of the mean 
level spacing) is only due to the finite ensemble used, 
whereas in reality the r.m.s. susceptibility diverges as the 
temperature becomes lower. To test this we compared 
the zero temperature CE r.m.s. susceptibility (which is 
identical to the EGCE susceptibility in this case) calcu- 
lated numerically on ensembles of different numbers of 
disorder realizations. We observed, indeed, that the nu- 
merically calculated r.m.s. susceptibility grows with the 
ensemble size, indicating it is really divergent, and thus 
temperature dependent for all temperatures. 

All the above refers to the CE or EGCE. Since the 
number of electrons is fixed, the probability of having a 
small separation between the lowest unoccupied level and 
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the highest occupied level is given by the Wigner distri- 
bution. In the GCE, however, there is a low probability 
for the fixed chemical potential to fall between two close 
level, and this probability is simply proportional to the 
separation itself when the separation is small. The tail of 
the susceptibility distribution thus varies as x -4 instead 
of x~ 3 m the CE or EGCE. This dependence is fitted to 
the GCE distributions in the insets of Fig. 03 and agrees 
quite well with the numerical results. Thus, the zero tem- 
perature r.m.s. susceptibility is expected to be finite in 
the GCE (although higher moments are still divergent). 
Numerical simulations have indeed shown that increasing 
the ensemble size only helps smoothing the r.m.s. GCE 
susceptibility. This means that the r.m.s. GCE suscep- 
tibility really saturates as a function of temperature for 
low enough temperatures. 

We thus arrive at a surprising conclusion: in contrast 
to what is usually taken for granted, the r.m.s. suscepti- 
bility is different in the CE (or EGCE) and in the GCE, 
diverging as a function of temperature when approaching 
absolute zero in the former case, while tending to a con- 
stant value in the former. Moreover, since the CE is qual- 
itatively similar to the EGCE at a lower temperature, the 
diverging temperature dependence explains why we can 
see a noticeable difference between the r.m.s. susceptibil- 
ities of the CE and the EGCE at very low temperatures. 
(It is not very large, however, since the temperature de- 
pendence is only logarithmic.) 

To conclude this section we may remark that since the 
main difference between the CE or EGCE and the GCE 
is connected with the events of small level separation at 
the Fermi level, we can understand why the average CE 
and EGCE susceptibilities show a temperature depen- 
dence even for temperatures below the mean level spac- 
ing. However, the average susceptibility is not divergent 
at zero temperature. This means that we cannot refer 
only to the tails of the susceptibility distribution, and 
the analysis becomes much more involved. In any case, 
however, we can conclude that the usual assumption of 
no temperature dependence when T < A cannot be taken 
for granted even when treating the average susceptibility. 



V. INTERACTION EFFECTS 

It is generally accepted that interaction should have a 
strong effect on the magnetic susceptibility. The inter- 
action contribution to the susceptibility is expected to 
be positive for a repulsive interaction and negative for 
an attractive oneS&i 9 *^. For diffusive systems this ef- 
fect is believed to be insensitive to the ensemble average 
chosen, and to have small fluctuations around its mean 
value. In these systems, for weak interaction and tem- 
peratures higher than the mean level spacing, this con- 
tribution varies as A|x L \(kp£) hxiTr/h) where r is the 
elastic mean free time and A is the dimensionless inter- 
action constant (The logarithmic factor is correct in the 
very low field regime $ < $ c = {T/E c ) 1/2 <S>o). For 



clean systems this contribution varies as A|x L \ (kpL) a for 
some geometry-dependent exponent a with a more com- 
plicated temperature-dependent factor. In the following, 
however, we will concentrate at zero temperature and 
repulsive interaction. 

In order to avoid the heavy computation involved in ex- 
act diagonalization or Hartree-Fock calculations, we will 
restrict ourselves to weak short-range interaction, which 
can be treated perturbatively to first order. We thus 
take as our model the usual Hubbard on-site interaction 
(electronic spin is restored in this section): 

U = C/ ff ^> s;T n s;i , (10) 

s 

where Uh specifies the interaction strength, and the ar- 
row subscripts denote spin direction. To first order in Uh 
the change in the system's energy at zero temperature is 
thus: 

AE = u H J2 n r( s ) n i(s), (ii) 

s 

where 

Tier 

n A s ) = ^2 \lpk;a(s)\ 2 
k=l 

is the number of electrons with spin projection a = f , j at 
lattice site s. In this later expression ipk-a- denotes the k- 
th single particle eigenfunctions without interaction, with 
spin projection a (The wave functions are, of course, spin 
independent; the spin index is used only for clarity); and 
n a is the total number of electrons with spin projection a 
in the system (which is conserved by the interaction). We 
will take these numbers as independent, as they may be 
set different experimentally using an in-plane magnetic 
field. This applies to the CE or EGCE; in the GCE the 
sum is over single-particle levels with energies lower than 
the Fermi energy (which can again depend on the spin 
direction) . 

The first-order interaction correction to the suscepti- 
bility is A% = —S~ 1 d 2 AE/dH 2 . Since for infinite clean 
systems n a (s) is position independent both with and 
without an applied magnetic, AE is field-independent 
and A% vanishes. This is in agreement with previous 
results—, according to which interaction corrections to 
the susceptibility are of the third order in the coupling 
strength (Uh in our case). On the other hand, for finite 
and/or disordered systems, n a (s) is non-uniform. More- 
over, an applied magnetic field breaks time-reversal sym- 
metry and drives the electronic density distribution to 
be more uniform. This makes AE smaller, and thus Ax 
is positive. This is again in accordance with previous 
resulta 7 i 8 i 19 . 20 . 

We now turn to numerical results. Since one can show 
that the change in the square of the absolute value of 
the wave-function amplitude at a given point due to the 
applied magnetic field is quadratic in the applied field, 
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FIG. 9: Corrections to the zero temperature susceptibility 
of spin-up electrons due to their interaction with spin-down 
electrons. The results are presented as a function of filling of 
spin-up electrons for different fillings of spin-down electrons 
in the CE (or EGCE), for both W = 2.0 (upper panel) and 
W = 4.0 (lower panel). 



we have 



Uh 
' S 



E 



d 2 ni 



dH 2 



-ni(s) +« T (s) 



dH 2 



(12) 



We may thus interpret the first term in the brackets as 
the correction to the orbital susceptibility of the spin-up 
electrons due to their interaction with spin down elec- 
trons, and similarly for the second term. The first term 
was calculated numerically with our usual parameters, 
and the results are shown in Fig. [21 for the two values of 
disorder strength (W = 2.0 and W = 4.0). Since we are 
using numerical differentiation, our results for W = 0, 
which seemed unreliable to us are not shown. For sim- 
plicity, We have taken the interaction strength as our 
unit of energy, i.e., Uh = t. 

As in the non-interacting case, for W = 2.0 we can 
see ballistic boundary- induced oscillations. The non- 



oscillating part of the interaction correction shows the ex- 
pected disorder dependence (varying as W~ 2 as I should) 
and is size independent (This latter result applies also 
to the oscillating part). An interesting feature is the 
strong peak in A^ when the filling of spin-up electrons 
equals the filling of spin-down electrons. This can be ex- 
plained by treating the interaction in the Hartree-Fock 
manner, as an effective single-particle potential. Obvi- 
ously, the interaction between two electrons in the same 
orbital single-particle level but different spins is stronger 
than in the case where the two electrons reside in two 
different orbital levels. Thus, the single-particle energies 
of spin-up electrons up to the n^-th level are raised by a 
repulsive interaction more than the single particle ener- 
gies above the n^-th level. Now, the expression for the 
van-Vleck susceptibility, Eq. |3 contains at zero temper- 
ature only energy denominators between occupied and 
unoccupied levels, and is thus dominated by the lowest 
unoccupied levels and the highest occupied levels. When 
the filling of spin-up electrons differs enough from the 
filling of spin-down electrons, the energies of these lowest 
unoccupied levels and highest occupied levels are raised 
by approximately the same amount by the interaction, so 
their differences do not change very much and the cor- 
rections to the susceptibility are small. However, when 
the filling of spin-up electrons is approximately equal to 
the filling of spin-down electrons, the highest unoccupied 
levels are raised by the interaction more than the lowest 
unoccupied levels so that the energy differences become 
smaller and is enhanced, as we have seen in the nu- 
merical results. 

To conclude this section we note that our numerical 
results indicate that the r.m.s. value of Ax is much 
larger than its average, and that the average of Ax in the 
GCE is significantly lower than its CE (or EGCE) value, 
in contrast to what is usually expected. It seems that 
the distribution of A\ at zero temperature shares with 
the distribution of the non-interacting susceptibility the 
property of having a long tail with diverging moments, 
causing a strong temperature and statistical-ensemble de- 
pendence. This point, however, deserves further investi- 
gation. 



VI. CONCLUSIONS 

To conclude this paper, we summarize our main find- 
ings : 

(i) For temperatures lower than the mean level spacing, 
we have shown how the requirement of gauge-invariance 
causes a strong cancellation between the Larmor and 
van-Vleck susceptibilities, and how it makes using usual 
averaging procedures, like RMT, inapplicable for our 
discussion. We have also seen that, contrary to what 
one usually expects, the r.m.s. susceptibility depends 
strongly both on temperature and statistical ensemble in 
this regime of low temperatures. This was explained in 
term of the long tail of the zero temperature suscepti- 
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bility distribution, resulting in a diverging r.m.s. zero 
temperature susceptibility in the CE or EGCE but not 
in the GCE, and causing a logarithmic temperature de- 
pendence in the former case for all temperatures down to 
zero. We believe this to be responsible also to the strong 
temperature dependence of the average CE and EGCE 
susceptibilities for this low temperature regime. 

(ii) For temperatures higher than the mean level spac- 
ing our numerical results generally agree with previous 
calculations. They, together with analytical calculations 
show that the difference between the CE and the EGCE, 
which was not considered so far in connection with or- 
bital susceptibility, is negligible for diffusive systems, and 
is not the leading order contribution in the clean limit. 
However, the meaning of the temperature where the 
EGCE susceptibility becomes equal to the GCE value, 
and even more - the temperature where the r.m.s. sus- 
ceptibility becomes equal to the average value, are not 
clear. Especially in the latter case, identifying the tem- 
perature with the Thouless energy is quite problematic. 

(iii) Considering interaction between the electrons, we 
have given a new interpretation of its first order con- 
tribution. We have also demonstrated its dependence on 
in-plane magnetic field, and explained this effect in terms 
of a single-particle picture. 
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APPENDIX: DIFFERENCE BETWEEN 
ENSEMBLES IN THE MESOSCOPIC REGIME 

In this appendix we apply the general formulation of 
Kamenev and Gefenii regarding the differences between 
ensembles to our problem. They calculated SF, the 
disorder-averaged difference between the canonical free 
energy F and the expression £l + fiN, the Legendre trans- 
formation of the grand canonical potential fl calculated 
for system-independent chemical potential fi. The num- 
ber of particles N is thus sample-specific in the grand- 
canonical case, and the chemical potential is set so that 
the disorder and thermal averaged particle number will 
be equal to the fixed particle number in the CE. Treating 
the regime of temperatures higher than the mean level 
spacing A, they give the following expression, correct to 
the zeroth order in the small parameter T ftp: 

ttT^ 1 /A\ w [°° t 2w ~ /A \ 

(A.l) 

where 

~ 1 f°° 

K{t) = — / dLoe- l ^/ A K(u) 

^ J — oo 



is the Fourier transform of the two-level correlation func- 
tion K{uj). The w — term is the difference between 
the EGCE and the GCE, while the w = 1 is the leading 
order difference between the CE and the EGCE. Only 
these terms will be considered in what follows. 

For temperatures lower than the Thouless energy Ec, 
the magnetic-field dependent part of K(w) can be ap- 
proximated by the zero-mode Cooperon contribution 
(since we are treating temperatures higher than the mean 
level spacing, perturbative analysis applies): 

A 2 ,2 _ , ,2 

where en ~ Ec($/&o) 2 is the leading order shift in the 
lowest Cooperon eigenmode due to the applied magnetic 
field (valid for small magnetic fields, whose flux through 
the sample $ is smaller than the quantum flux <f>o). 

Inserting this into Eq. (| A. 1|> we get for the difference 
between the EGCE and the GCE the following expres- 
sion, previously derived by Oh et al&, for the magneti- 
zation per unit area: 

AW - A/gce ~ ^ f dt— * e-^m*. 

1 ®o J-oo smh (nt) 

(A.3) 

From this it follows that 

A/egce - A/gce ~ j=- r * *° ' Z * " (A.4) 

where $ c = (T/E c ) 1 <J> . The difference between the 
CE and the EGCE is given by: 

(A. 

From this it follows that 




In our parameter regime we can see the the difference 
between the CE and the EGCE is parametrically smaller 
than the difference between the EGCE and the GCE. In 
fact, the former effect is of the same order of magnitude as 
the contribution of higher order perturbative corrections 
to K(w) to the latter. This is in accordance with our 
numerical results. 
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